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The charged lepton-antilepton pair signature, in par- 
ticular e + e~ and /i + ^t~, has been a leading discovery 
channel for new particles such as the J/tp and T mesons 
and the Z boson since they have cleaner experimental 
signatures and lower backgrounds than hadronic signa- 
tures. 

Many models beyond the standard model (SM) pre- 
dict the existence of new particles decaying to lepton- 
antilepton pairs. The Eq Z's [l| and the Randall- 
Sundrum (RS) graviton are examples of specific new 
particles decaying to a lepton-antilepton final state. The 
Z^p, Z' x , Z' v , Zj, Z' sec , and Z' N are chosen for testing the 
Eq model. Assuming one extra dimension, we test the 
RS model in the k/Mpi range between 0.01 and 0.1 [H , 
where k is the curvature of the extra dimension and M pi 
is the reduced effective Planck scale. 

In recent publications, the CDF [3] and DO [H] Collab- 
orations set limits on these models with 1.3 and 1.0 fb" 1 
of integrated luminosity and limits on Z' with SM cou- 
pling and RS graviton for fc/Mp;=0.1 are 923 and 900 
GeV/c 2 , respectively. Using a data set twice as large (2.5 
fb _1 ), this Letter describes a search for e + e~ resonances 
in the invariant mass range of 150-1,000 GeV/c 2 , and 
we set upper limits on <r(pp — > X) ■ B(X — > e + e _ ) at 
the 95% credibility level (C.L.) where X is a spin 1 or 
spin 2 particle. We also set lower mass bounds on the Z' 
with SM coupling, the Z's in the E§ model, and the RS 
graviton. 

This analysis is based on data collected with the CDF 
II detector [6j. The relevant components of the detec- 
tor for this analysis are the tracking system and the 
calorimeters. The tracking system consists of a 96 layer 
drift chamber called the central outer tracker (COT), 
surrounding an eight-layer silicon tracker. Both are in- 
side a 1.4 T solenoidal magnet. The COT covers the 
range of pseudorapidity |ry| <1.1 0, and the silicon 
tracker covers |?7| up to 2.0. The electromagnetic (EM) 
and hadronic calorimeters, which are sandwiches of lead 
(EM) or iron (hadronic) absorber and plastic scintilla- 
tor. They are outside the magnet, and are divided into 
a central calorimeter (|?7| <1.1) and two plug calorime- 
ters (l.K \rj\ <3.6). Both the central and the plug EM 
calorimeters have fine-grained shower profile detectors at 
EM shower maximum. 
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We use the same on-line event selection criteria (trig- 
gers) used in our previous report Off-line events are 
required to have two isolated electrons Q, one in the 
central EM calorimeter and the other one in either the 
central (CC) or the plug (CP) EM calorimeters. Only 
electrons with Et Q greater than 25 GeV and M < 2 are 
used in order to ensure 100% trigger efficiency and cover- 
age by the the silicon tracker. Electrons in the central EM 
calorimeter are required to have a well-measured track 
in the COT system pointing at an energy deposit in the 
calorimeter. For electrons in the plug EM calorimeter, 
the track association uses a calorimeter-seeded silicon- 
tracking algorithm An opposite-charge requirement 
is applied to electron-objects pairs detected in the central 
EM calorimeter. No such requirement is applied when 
one electron is detected in the plug, where ^-dependent 
charge misidentification occurs. Events with both elec- 
trons in the plug EM calorimeter are not considered in 
this Letter since adding them gains little sensitivity. 

The pythia [l(| Tune A [ll| Monte Carlo event gen- 
erator is used to model the expected signals and back- 
grounds unless otherwise stated. For spin 1 Z' , SM-like 
couplings are assumed, and for spin 2 resonances, the 
RS graviton model with k/Mpi=0.l is used. The total 
selection efficiencies of spin 1 particles vary from 27 to 
38% and those of spin 2 particles vary from 28 to 32% as 
functions of the particle mass in the search range. 

There are three sources of background. One is Drcll- 
Yan production of e + e _ pairs (DY), which is the dom- 
inant source of background and is irreducible. Another 
is dijets and VF+jets production (referred to as "QCD" 
background) where one or more jets is misidentified as 
electron. Other contributions include Z/j* — > t + t~, 
ti, and diboson (W7, WW, WZ, ZZ, 77) production that 
collectively are referred to as "other SM" backgrounds. 

The simulated DY prediction is normalized to the data 
after subtracting other SM and QCD backgrounds in an 
invariant mass window from 76 to 106 GeV/c 2 for CC 
events and from 81 to 101 GeV/c 2 for CP events to es- 
timate the DY background. Different mass windows are 
used because the QCD background rate in CP events is 
higher than in CC events. We assign a 3.6% systematic 
uncertainty in the DY prediction to take into account the 
invariant-mass dependence of the fc-factor [13[ that is the 
difference between the leading and the next-to-next-to- 
leading order DY cross sections. The uncertainty in the 
DY prediction due to the choice of the parton distribution 
function set CTEQ6M [3 using the Hessian method OH 
is 3.7-6.4-13% (200-600-1,000 GeV/c 2 ) depending on 
the invariant mass. 

The QCD background estimation is determined from 
the experimental data. The estimate is obtained us- 
ing the probability for a jet to be misidentified as an 
electron [3]. We measure this probability with a jet- 
triggered data sample. We then apply the misidentifi- 
cation probability to each jet in events with one good 
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electron candidate and one or more jets. To estimate the 
dijet background contribution, events with W or Z can- 
didates are removed from the sample before applying the 
jet misidentification probability (MP). Events with W 
candidates are identified with one goo d electron and a 
large missing transverse energy $t [171 ] and events with 
Z candidates are identified with two "loose electrons". 
To estimate the VF+jets background, events with Z can- 
didates are removed and events with W candidates are 
retained. The dominant systematic uncertainty in the 
predicted QCD background is due to the 20% uncertainty 
in the jet MP, which is obtained from the variation in the 
MP measured in the different jet data samples. 

Other SM contributions to the background are esti- 
mated with simulation samples [12J. These simulated 



samples are normalized to the product of the theoreti- 
cal cross sections and the integrated luminosity. Fig. Q] 
shows the observed e + e~ invariant mass spectrum from 
2.5 fb _1 of data together with the expected backgrounds. 



affects the shape of the e + e~ invariant mass distribu- 
tion. The uncertainty on the acceptance due to parton- 
distribution-function uncertainties is evaluated using the 
same method that was used for the DY prediction, and 
found to be 1.9% for CC and 0.6% for CP events. 

The search for e + e~ resonances in the high-mass range 
of 150-1,000 GeV/c 2 uses an unbinned likelihood ratio 
statistic, A, defined in Eqs.[T]-[5] 20]: 



max. Lb 

n b >0 



max. C s +b 

n b >0,n s >0 



< A < 1, < -2 In A < oo (1) 



(n s + n b ) N e-^+ n ^ A n,S(xi\n) + n b B( Xl ) 

L s+b ~ —, 1 1 — 



n s + rib 



N 



N\ 



Y[B( Xi ). 



(2) 
(3) 
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FIG. 1: Invariant mass distribution of e + e~ events compared 
to the expected backgrounds. Dots with error bars are data. 
The dark shaded region represents "other SM" background, 
the light shaded region shows "QCD" background, and the 
white region corresponds to Zj^j* — > e + e~ background. The 
inset shows the same for the 240 GeV/c 2 region. The hatched 
histogram shows the shape of the expected signal from a 240 
GeV/c 2 spin 1 particle (of negligible intrinsic width) on top 
of the total background. The hatched region is normalized to 
the number of excess events seen in the data. 



where Lb is the likelihood for a null hypothesis that is de- 
scribed by the SM only, C s +b is the likelihood for a test 
hypothesis that is described by physics beyond the SM 
together with the SM. The quantities n s and n b are the 
number of signal and background candidates which are 
determined by the fit and N is the number of candidates 
observed in data, each represented by a vector {xi} of ob- 
servables. The signal probability density function (PDF), 
S(x\fi), is a Gaussian with a floating mean [i and a fixed 
width, and B(x) is a background PDF obtained from 
the total background template. The widths of the sig- 
nal PDF are determined from simulation {(Jm„ — 0.8565 
GeV/c 2 +0.0192 • M ee for M ee > 150 GeV/c 2 ) with the 
assumption that the decay widths of resonances are much 
less than the experimental resolution. The quantities 
C s +b and the Cb are maximized separately without exter- 
nal background constraints. The function —2 In A is cal- 
culated over the search range of 150-1,000 GeV/c 2 and 
the most prominent local maxima are listed in Table [J 
The most significant deviation between data and the SM 
prediction occurs at an invariant mass of 241.3 GeV/c 2 



where — 21nA is 14.4. The (data — background)/^ [21 ] 
corresponding to the region of maximum —2 In A is cal- 
culated by counting the number of observed events and 
estimated backgrounds within ±2 (Jm„ of the maximum, 
and it is 3.8. 



The systematic uncertainty for other SM backgrounds 
is dominated by the 6% uncertainty in the integrated 
luminosity measurement 11811 and 8% uncertainty in the 
theoretical cross sections [19|. Other systematic sources 
are the uncertainty on the scale factor of electron identifi- 
cation efficiency that comes from the difference between 
data and simulated events (1.3% for CC and 2.3% for 
CP events), the energy scale (1.0%), and the energy res- 
olution (0.6% for CC and 0.3% for CP events), which 



TABLE I: The prominent local maxima in the search range 
of 150-1,000 GeV/c 2 . 

M x (GeV/c^) 241.3 272.7 478.9 725.2 
— 2 In A 14.4 3.7 2.6 4.1 



To estimate the probability of observing an excess 
equal to or greater than the maximum observed ex- 
cess anywhere in the search range of 150-1,000 GeV/c 2 , 
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we simulated 100,000 experiments assuming background 
only. The distribution of maximum —2 In A on these sim- 
ulated experiments is shown in Fig.O Assuming only SM 
physics, the probability of observing a number of events 
equal to or greater than the observed excess is defined 
as the fraction of simulated experiments with maximum 
—2 In A equal to or greater than 14.4, and is 0.6% which 
corresponds to the 2.5 a level of excess over the back- 
ground. 



Max.-2liU=14.4 




20 25 30 

Maximum -2lnA. 



FIG. 2: Distribution of maximum —2 In A in simulated experi- 
ments that assume only background. The arrow indicates the 
value observed in data: — 21nA=14.4. 

Upper limits on a(pp — > X) ■ B(X — > e + e~) at the 
95% C.L. are calculated as a function of mass using a 
Bayesian binned likelihood method with a full consider- 
ation of uncertainties on the total signal efficiency and 
the background estimation (22|. Fig. [3] (a) shows the ob- 
served upper limits from data and the expected limits 
from background-only simulated events for spin 1 parti- 
cles as a function of the e + e~ invariant mass, together 
with the expected cross sections for Z's Fig. [3] (b) 
shows the same but for spin 2 particles, together with the 
expected cross sections for RS gravitons. The cross sec- 
tions for Z's and RS gravitons are calculated at leading 
order with pythia and then multiplied by a factor of 1.3 
in order to approximate a next-to-leading-order predic- 
tion as done in reports of earlier results. Table ITT1 shows 
the lower mass limits of the SM coupling and E§ Z's and 
Fig. [4] shows the excluded RS graviton mass region with 
respect to k/Mpi. 



TABLE II: Expected and observed 95% C.L. lower limits on 
Z's masses. 



Z' Model 


Z'sM 


z'^ 


Z' x Z' n Z'j 


7' 

^ sec 


Z'n 


Expected Limit (GeV/c 2 ) 
Observed Limit (GeV/c 2 ) 


961 


846 


857 873 755 


788 


831 


963 


851 


862 877 735 


792 


837 



To conclude, we have searched for e + e resonances 
with 2.5 fb _1 of data collected by the CDF II detector. 
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FIG. 3: The upper limits on a(pp -> X) ■ B(X -» e + e") as 
function of the mass of an X particle at the 95% C.L. where 
X is a spin 1 particle (a) or a spin 2 particle (b) together with 
model predictions. 



The largest excess over the standard model prediction is 
at an e + e~ invariant mass of 240 GeV/c 2 . The probabil- 
ity of observing such an excess arising from fluctuation 
in the standard model anywhere in the mass range of 
150-1,000 GeV/c 2 is 0.6%. We also set upper limits on 
a(pp -> X) ■ B{X -> e+e") at the 95% C.L. for spin 
1 and spin 2 particles. The SM coupling Z' with mass 
below 963 GeV/c 2 and the Eg Z's with masses below 
735/877 (lightest/heaviest) GeV/c 2 are excluded at the 
95% C.L. RS gravitons with masses below 848 GeV/c 2 
are excluded at the 95% C.L. for k/M P i =0.1. 
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